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Accessible quantification of multiparticle entanglement 
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Entanglement is a key ingredient for quantum technologies and a fundamental signature of quantumness 
in a broad range of phenomena encompassing many-body physics, thermodynamics, cosmology, and life sci¬ 
ences. For arbitrary multiparticle systems, entanglement quantification typically involves nontrivial optimisation 
problems, and may require demanding tomographical techniques. Here we develop an experimentally feasible 
approach to the evaluation of geometric measures of multiparticle entanglement. Our framework provides ana¬ 
lytical results for particular classes of mixed states of N qubits, and computable lower bounds to global, partial, 
or genuine multiparticle entanglement of any general state. For global and partial entanglement, useful bounds 
are obtained with minimum effort, requiring local measurements in just three settings for any N. For genuine 
entanglement, a number of measurements scaling linearly with N is required. We demonstrate the power of our 
approach to estimate and quantify different types of multiparticle entanglement in a variety of Wqubit states 
useful for quantum information processing and recently engineered in laboratories with quantum optics and 
trapped ion setups. 


I. INTRODUCTION 

The fascination with quantum entanglement has evolved 
over the last eight decades, from the realm of philosophi¬ 
cal debate HI to a very concrete recognition of its resource 
role in a range of applied sciences ElEl. While considerable 
progress has been achieved in the detection of entanglement 
HHia, its experimentally accessible quantification remains 
an open problem for any real implementation of an entangled 
system OMl. Quantifying entanglement is yet necessary to 
gauge precisely the quantum enhancement in information pro¬ 
cessing and computation EiiaiMi, and to pin down exactly 
how much a physical or biological system under observation 
departs from an essentially classical behaviour 1^ . This is 
especially relevant in the case of complex, multiparticle sys¬ 
tems, for which only quite recently have notable advances 
been reported on the control of entanglement Il26i429ll . 

An intuitive framework for quantifying the degree of mul¬ 
tiparticle entanglement relies on a geometric perspective Il30l - 
\32\ . Within this approach, one first identifies a hierarchy 
of non-entangled multiparticle states, also referred to as M- 
separable states for 2 < M < N, where N is the number of 
particles composing the quantum system of interest; see Fig- 
ure[2 Introducing then a distance functional D respecting nat¬ 
ural properties of contractivity under quantum operations and 
joint convexity (see Methods) m, the quantity defined 
as 

Em(q)^ inf D(e,g), (1) 

q M-separable 

is a valid geometric measure of (M-inseparable) multiparticle 
entanglement in the state g. 

Some special cases are prominent in this hierarchy. For 
M - N, the distance from A-separable (also known as fully 


* cianciamso.marco@gmail.com 
’ thomas.r.bromley@gmail.com 
^ gerardo.adesso@nottingham.ac.uk 



FIG. 1. Geometric picture of multiparticle entanglement in a quan¬ 
tum system of N particles. Each red-shaded convex set contains M- 
separable states (for 2 < M < N), which can be defined as follows. 
The set of pure M-separable states is given by the union of the sets of 
tensor products |Tm) = li/'®) of pure states |i^®), with respect 

to any partition of the N particles into M subsystems k = 1,..., M; 
the set of general (mixed) M-separable states is then formed by all 
convex mixtures of pure M-separable states, where each teim in the 
mixture may be factorised with respect to a different multipartition. 
This provides a partition-independent classification of separability 
(see also the Supplementary Material). For any M, the multiparti¬ 
cle entanglement measure of a state g is defined as the minimum 
distance, with respect to a contractive and jointly convex distance 
functional D, from the set of M-separable states. We refer to the case 
M = N (dotted line) as global entanglement, the case M = 2 (solid 
line) as genuine entanglement, and any intermediate case (dashed 
line) as partial entanglement, as detailed in the main text. 


separable) states defines the global multiparticle entanglement 
which accounts for any form of entanglement distributed 
among two or more of the N particles. Geometric measures 
of global entanglement have been successfully employed to 
characterise quantum phase transitions in many-body systems 
GS and directly assess the usefulness of initial states for 






2 


Grover’s search algoritms ll35l . On the other extreme of the 
hierarchy, for M - 2, the distance from 2-separable (also 
known as biseparable) states defines instead the genuine mul¬ 
tiparticle entanglement E!^, which quantifies the entanglement 
shared by all the N particles, that is the highest degree of 
inseparability. Genuine multiparticle entanglement is an es¬ 
sential ingredient for quantum technologies including mul¬ 
tiuser quantum cryptography ||36|, quantum metrology llJTl . 
and measurement-based quantum computation 138]. Finally, 
for any intermediate M, we can refer to as partial multi¬ 
particle entanglement. The presence of partial entanglement 
is relevant in quantum informational tasks such as quantum 
secret sharing Col and may play a relevant role in biological 
phenomena Il25ll39]l . Probing and quantifying different types 
of entanglement can shed light on which nonclassical fea¬ 
tures of a mixed multiparticle state are necessary for quantum- 
enhanced performance in specific tasks 13 and can guide the 
understanding of the emergence of classicality in multiparticle 
quantum systems of increasing complexity go). 

The quantitative amount of multiparticle entanglement, be 
it global or genuine (or any intermediate type), has an intu¬ 
itive operational meaning when adopting the geometric ap¬ 
proach. Namely, E’^ measures how distinguishable a given 
state Q is from the closest M-separable state. Given some 
widely adopted metrics, such a distinguishability is directly 
connected to the usefulness of g for quantum information pro¬ 
tocols relying on multiparticle entanglement. For instance, the 
trace distance of entanglement is operationally related to the 
minimum probability of error in the discrimination between 
Q and any M-separable state with a single measurement 133]. 
Furthermore, the geometric entanglement with respect to rela¬ 
tive entropy or Bures distance sets quantitative bounds on the 
number of orthogonal states that can be discriminated by local 
operations and classical communication (LOCC) HTl . The 
geometric entanglement based on infidelity ll3TI (monotoni- 
cally related to Bures distance) has also a dual interpretation 
based on the convex roof construction Ba, that is, it quanti¬ 
fies the minimum price (in units of pure-state entanglement) 
that has to be spent on average to create a given density matrix 
p as a statistical mixture of pure states. 

It is therefore clear that finding the minimum in Eq. 
and hence evaluating geometric measures of multiparticle en¬ 
tanglement defined by meaningful distances, is a central chal¬ 
lenge to benchmark quantum technologies. However, obtain¬ 
ing such a solution for general multiparticle states is in prin¬ 
ciple a formidable problem. Even if possible, there would 
remain major challenges for experimental evaluation, which 
would in general require a complete reconstruction of the state 
through full tomography. Eor multiparticle states of any rea¬ 
sonable number of qubits, full state tomography places signif¬ 
icant demands on experimental resources, and it is thus highly 
desirable to provide quantitative guarantees on the geometric 
multiparticle entanglement present in a state, via non-trivial 
lower bounds, in an experimentally accessible way lfT3l423l . 

Here we provide substantial advances towards addressing 
this problem in a general fashion. We identify a general 
framework for the provision of experimentally friendly quan¬ 
titative guarantees on the geometric multiparticle entangle¬ 


ment present in a state. This approach consists of: 

(1) Choosing a set of reference states: Eind a restricted fam¬ 

ily of N-qubit states with the property that any state may 
be mapped into this family through a fixed procedure 
of single-qubit LOCC. This reference family should be 
simple to characterise, and can be chosen from experi¬ 
mental or theoretical considerations. 

(2) Identifying M-separable reference states: Apply the 

fixed LOCC procedure to the general set of M- 
separable states, hence identifying the subset of 
M-separable states within the reference family. 

(3) Calculating E'^ for the reference states: Solve the opti¬ 

misation problem for the geometric entanglement of 
reference states. This is dramatically simplified by us¬ 
ing the properties of contractivity and joint convexity, 
that hold for any distance functional D defining a valid 
entanglement measure, and imply in particular that one 
of the closest M-separable states to any reference state 
is to be found itself within the reference family. 

(4) Deriving optimised lower bounds for any state: Exploit 

the freedom to apply single-qubit unitaries to any 
A-qubit state g in order to find the corresponding ref¬ 
erence state with the highest geometric entanglement, 
providing an optimised lower bound to E^{g). 

This process presents a versatile and comprehensive approach 
to obtain lower bounds on geometric multiparticle entangle¬ 
ment measures according to any valid distance. While build¬ 
ing on some previously utilised methods for steps (1) GMa 
and (4) lfT3l fBl [T^. it introduces novel techniques in steps 
(2) and most importantly (3), which are crucial for completing 
the framework and making it effective in practice (see e.g. Ap¬ 
pendix C, D, E, and E of the Supplementary Material). 

To illustrate the power of our approach, we focus initially 
on a reference family of mixed states nr of A qubits, that we 
label states, which form a subset of the class of states 
having all maximally mixed marginals. This family includes 
maximally entangled Bell states of two qubits and their mix¬ 
tures, as well as multiparticle bound entangled states Il43lj46l . 
For any A, these states are completely specified by three eas¬ 
ily measurable quantities, given by the correlation functions 
Cj - (cr®^), where {crj}j^i 2 ,i are the Pauli matrices. In the 
following we show how every entanglement monotone E^ 
can be evaluated exactly for any even A on these states, by 
revealing an intuitive geometric picture common to all valid 
distances D. For odd A, the results are distance-dependent; 
we show nonetheless that can still be evaluated exactly 
if D denotes the trace distance. The results are nontrivial for 
all M > \NI2'\ in the hierarchy of Figure [T] A central ob¬ 
servation, in line with the general framework, is that an arbi¬ 
trary state of A qubits can be transformed into an state by 
a LOCC procedure, which cannot increase entanglement by 
definition. This implies that our exact formulae readily pro¬ 
vide practical lower bounds to the degree of global and partial 
multiparticle entanglement in completely general states. Im¬ 
portantly, the bounds are obtained by measuring only the three 
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correlation functions {c/) for any number of qubits, and can 
be further improved by adjusting the local measurement basis 
(see Figure]^ for an illustration). 

Furthermore, we discuss how our results can be extended 
to allow for the quantitative estimation of genuine multipar¬ 
ticle entanglement as well, at the cost of performing extra 
measurements. Since states are always biseparable, we 
must consider a different reference family. We focus on the 
class of A^-qubit states obtained as mixtures of Greenberger- 
Horne-Zeilinger (GHZ) states ElSii, the latter being central 
resources for quantum communication and estimation; this 
class of states depends on 2^ - 1 real parameters. We cal¬ 
culate exactly distance-based measures of genuine multiparti¬ 
cle entanglement for these states, for every valid D. Once 
more, these analytical results provide lower bounds to geo¬ 
metric measures of genuine entanglement for any general state 
of N qubits, obtainable experimentally in this case by per¬ 
forming at least H- 1 local measurements EH. 

We demonstrate that our results provide overall accessible 
quantitative assessments of global, partial, and genuine multi¬ 
particle entanglement in a variety of noisy states produced in 
recent experiments Il40ll44l 149145^ . going beyond mere detec¬ 
tion EHI2, yet with a signficantly reduced experimental over¬ 
head. Compared with some recent complementary approaches 
to the quantification of multiparticle entanglement 1(13142^ . 
we find that our results, obtained via the general quantitative 
framework discussed above, fare surprisingly well in their ef¬ 
ficiency and versatility despite the minimal experimental re¬ 
quirements (see Table|^for an in-depth comparison). 


II. RESULTS 

A. Global and partial multiparticle entanglement 

We begin by choosing as our reference family the set 
of A-qubit A4^ states. An state is defined as vj - 
^ CjCrJ^^, where I is the 2x2 identity matrix. 

These states are invariant under permutations of any pair of 
qubits and enjoy a nice geometrical representation in the space 
of the three correlation functions cj, corresponding to a tetra¬ 
hedron for even N and to the unit ball for odd N, as depicted in 
Figure]^ We can then characterise the subset of M-separable 
Ai^ states for any N. We find that, if M > rA^/2], then the 
M-separable Ad^ states hll a subset corresponding to an octa¬ 
hedron in the space of the correlation functions (see Figure]^. 
When M < \NI2\ all Ad^ states are instead M-separable. The 
proofs are deferred to the Supplementary Material. 

We can now tackle the quantification of global and partial 
multiparticle entanglement in these states; for the latter, we 
will always focus on the nontrivial case M > rA^/21 through¬ 
out this section. First, we observe that the closest M-separable 
state to an Ad^ state m, which solves the optimisation in 
Eq. Q, can always be found within the subset of M-separable 
Ad^ states, yielding a considerable simplihcation of the gen¬ 
eral problem. To hnd the exact form of and consequently 
of we approach the cases of even and odd N sepa¬ 


rately. 

For even N, we prove that there exists a unique solution 
to the minimisation in Eq. independent of the specific 
choice of contractive and jointly convex distance D. Namely, 
the closest M-separable state g^ is on the face of the octa¬ 
hedron bounding the corner of the tetrahedron in which in¬ 
is located, and is identihed by the intersection of such octa¬ 
hedron face with the line connecting nr to the vertex of the 
tetrahedron corner, as depicted in Figure [^d). It follows that, 
for any nontrivial M, valid D, and even N, the multiparticle 
entanglement £’^(iir*‘'r*) of an Ad^ state with correlation 
functions (cj) is only a monotonically increasing function of 
the Euclidean distance between the point of coordinates (cj) 
and the closest octahedron face, which is in turn proportional 
to hj^ - \cj\ — 1) (notice that equals the bipartite 

measure known as concurrence for N - 2 1231531). We 
have then a closed formula for any valid geometric measure 
of global and partial multiparticle entanglement on an arbi¬ 
trary Ad^ state with even N, given by 




0, < 0 (or M < A/2); 

fnihu), otherwise. 


( 2 ) 


where fo denotes a monotonically increasing function whose 
explicit form is specific to each distance D. In Table [I^ we 
present the expression of fo for relevant distances in quantum 
information theory. 

Eor odd N, the closest M-separable state ^^7 to any Ad^ 
state nr is still independent of (any nontrivial) M. However, 
different choices of D in Eq. o are minimised by different 
states g^. We focus on the important but notoriously hard-to- 
evaluate case of the trace distance Djfm, g) (see Table [Till. In 
the representation of Eigure [^c), the trace distance amounts 
to half the Euclidean distance on the unit ball. It follows that 
the closest M-separable state to m is the Euclidean or¬ 
thogonal projection onto the boundary of the octahedron, see 
Eigure |^e). We can then get a closed formula for the trace 
distance measure of global and partial multiparticle entangle¬ 
ment of an arbitrary Ad^ state with odd N as 

well, given by 




0, Act < 0 (or M < rAt/2]); 
^ 0 <Act<3|c,-|/2V/; 

min -H 5 ( 2 Act - \cj\)^ , 


otherwise. 


(3) 

The usefulness of the just derived analytical results for mul¬ 
tiparticle entanglement is not limited to the Ad^ states. In ac¬ 
cordance with our general framework, a crucial observation is 
that the Ad^ states are extremal among all quantum states with 
given correlation functions (c^). Specifically, any general state 
Q of N qubits can be transformed into an Ad^ state with the 
same {cj) by means of a procedure that we name Ad^-fication, 
involving only LOCC (see Methods). This immediately im¬ 
plies that, for any M > \NI2\ the multiparticle entanglement 
of Q can have a nontrivial exact lower bound given by the 
corresponding multiparticle entanglement of the Ad^ state m 
with the same (c/). 


£S(u 7'^^>) < £S(p), Vp : Tr(pa-®^) - (/ = 1,2,3). (4) 
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Comparison Table: Experimentally friendly methods to quantify multiparticle entanglement 


Ref. 

Experimental 

friendliness 

Computational 

friendliness 

M-inseparability 

quantified 

Entanglement 

measure 

caiaa 

Variable 

Optimisation required 

2 < M < A 

Any convex and continuous measure 

(m 

0(2^) 

Optimisation required 

2 

Genuine multiparticle concurrence 

CHIIol 

0(A) 

Closed formula 

2 

Robustness of entanglement 

(H 

0(A) 

Closed formula 

2 

Genuine multiparticle concurrence 

(m 

0(A) 

Closed formula 

2 

Genuine multiparticle concurrence 

[I5l[l6l 

0(A) 

Closed formula 

2 

Polynomial invariant (three-tangle) 

ca 

0(A) 

Closed formula 

2 

Genuine multiparticle negativity 

03 

0(A) 

Closed formula 

2 < M < A 

Infidelity-based geometric measure 

[*] (GHZ-diagonal set) 

0(A) 

Closed formula 

2 

All distance-based measures 

[*] {Ml set) 

3 

Closed formula 

rA/21 < M < A 

All distance-based measures 


TABLE I. A comparison of relevant literature on experimentally friendly quantification of multiparticle entanglement (based on accessible 
lower bounds). For each reference, we give the experimental friendliness, in terms of the number of local measurement settings required, and 
also the computational friendliness. The levels of M-inseparability quantified are given, along with the entanglement measures to which each 
work applies; [>i<] refers to this paper. 


Distance D 


foihr.) 

Relative 
entropy Dre 

Tr[£)(log2£)-log2^)] 

i[(l-Mlog2(l-M 

+ (1 +hj^)\og^{\ +h^)\ 

Trace Dj, 

^Trlp-^l 


Infidelity Dp 

l-[Tr(^V?£>V?)J 

f(i-VT^) 

Squared 
Bures Dp 

2|l-Tr(VV?eV?)| 

2 - VT^- vm; 

Squared 
Hellinger Dh 

2[l-Tr(VeV?)] 

2 - vm; 


TABLE 11. Analytical expression of global and partial multiparticle 
entanglement for states of an even number N of qubits as 
defined by Eq. |^, for representative choices of the distance D. 

From a practical point of view, one needs only to measure the 
three correlation functions (Cy), as routinely done in optical, 
atomic, and spin systems naiaiigEoi, to obtain an estimate 
of the global and partial multiparticle entanglement content of 
an unknown state g with no need for a full state reconstruction. 

Furthermore, the lower bound can be improved if a partial 
knowledge of the state g is assumed, as is usually the case for 
experiments aiming to produce specific families of states for 
applications in quantum information processing IBlIlsiiig. 
In those realisations, one typically aims to detect entangle¬ 
ment by constructing optimised entanglement witnesses tai¬ 
lored on the target states ||4l. By exploiting similar ideas, we 
can optimise the quantitative lower bound in Eq. 0 over all 
possible single-qubit local unitaries applied to the state g be¬ 
fore the Af^-fication, 

sup£'^(?ir'‘^^*) < E^{Ug,gUl) = E^{g ), (5) 

t/® 

where Tr(t/,8pt/g cr®^) - cj and = ^a=i denotes a 
single-qubit local unitary operation. Experimentally, the opti¬ 
mised bound can then be still accessed by measuring a triple 
of correlations functions {cj} given by the expectation val¬ 


ues of correspondingly rotated Pauli operators on each qubit, 
Cj - as illustrated in Eigurej^a). Optimality 

in Eq. ([^ can be achieved by the choice of t/® such that 
- 5 ( 2 |C;I - 1) is maximum. The optimisation proce¬ 
dure can be significantly simplified when considering a state 
g which is invariant under permutations of any pair of qubits. 
In such a case, one may need to optimise only over three an¬ 
gles {fl, if/, (f)} parameterising a generic unitary applied to each 
single qubit; the optimisation can be equivalently performed 
over an orthogonal matrix acting on the Bloch vector of each 
qubit (see Methods). 

We can now investigate how useful our results are on con¬ 
crete examples. Table[III|presents a compendium of optimised 
analytical lower bounds on the global and partial multiparti¬ 
cle entanglement of several relevant families of A-qubit states 
Il46ll47llfni52ll54ti58ll62ti64l. up to A = 8. All the bounds 
are experimentally accessible by measuring the three correla¬ 
tion functions (c;), corresponding to optimally rotated Pauli 
operators (see also Eigure]^. 

Let us comment on some cases where our analysis is par¬ 
ticularly effective. Eor GHZ states, cluster states, and half- 
excited Dicke states, which constitute primary resources for 
quantum computation and metrology El [381, we get the 
maximum hj^ - \ for any even A. This means that our bounds 
remain robust to estimate global and partial entanglement in 
noisy versions of these states (i.e. when one considers mix¬ 
tures of any of these states with probability q and the max¬ 
imally mixed state with probability I - q) for all q > 1/3. 
Notably, for values of q sufficiently close to 1, our bounds 
to global entanglement can be tighter than the (more experi¬ 
mentally demanding) ones derived very recently in Ref lfT3l . 
as shown in Eigure|^a). Eocusing on noisy GHZ states, we 
observe however that our scale-invariant threshold q > 1/3, 
obtained by measuring the three canonical Pauli operators for 
each qubit, is weaker than the well-established inseparabil¬ 
ity threshold q > 1/(1 -t 2^^') S. Nevertheless, we note 
that our simple quantitative bound given by Eq. 0 becomes 
tight in the paradigmatic limit of pure GHZ states {q - 1) of 
any even number A of qubits, thus returning the exact value 
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Alice Bob Charlie Natalie 



FIG. 2. Experimentally friendly protocol to quantify global and partial At-particle entanglement. Top row: (a) A state goiN qubits is shared by 
N parties, named Alice, Bob, Charlie, ..., Natalie. Each party, labelled by a = A,..., A, locally measures her or his qubit in three orthogonal 
directions with j = 1,2,3, indicated by the solid arrows. If the shared state g is completely unknown, a standard choice can be to 

measure the three canonical Pauli operators for all the qubits (corresponding to the directions of the dashed axes); if instead some partial 
information on g is available, the measurement directions can be optimised a priori. Once all the data are collected, the N parties communicate 
classically to construct the three correlation functions {cy), with Cj = d"*"*). Middle row: For any A, one can define a reference subset of 

A-qubit states with all maximally mixed marginals (Al^ states), which are completely specified by a triple of orthogonal correlation functions 
{Cjj. These states enjoy a convenient representation in the space of jci,C 2 ,C 3 ). (b) For even A, AlJ, states fill the tetrahedron with vertices 
{1, (-1)^^^, 1), {-1, -(-1)"/^, 1), {1, -(-1)'''^^, -1) and {-1, (-1)"/^, -1). (c) For odd A, they are instead contained in the unit Bloch ball. For 
any M > rA/2], M-separable Af^ states are confined to the octahedron with vertices j±l,0,0), {0, ±1,0) and {0,0, ±1), illustrated in red in 
both panels; conversely, for M < rA/2], all Af^ states are M-separable. Bottom row: Geometric analysis of multiparticle entanglement. The 
bottom panels depict zooms of (d) a corner of the tetrahedron for even A and (e) a sector of the unit sphere for odd A, opposing a face of the 
octahedron of M-separable Al^ states (for M > rA/2]). Instances of inseparable states are indicated by blue circles, and their closest 
M-separable states by red crosses. The cyan surfaces in each of the two bottom panels contain states with equal global and partial multiparticle 
entanglement which we compute exactly. The results are valid for any contractive and jointly convex distance D in the even A case, and 
for the trace distance in the odd A case. The entanglement of an Al^ state with correlation functions {ci,C 2 ,C 3 l provides an analytical lower 
bound for the entanglement of any A-qubit state with the same correlation functions, such as the state g initially shared by the A parties in (a). 
The bound is effective for the most relevant families of A-qubit states in theoretical and experimental investigations of quantum information 
processing, as we show in this article. 


of their global multiparticle entanglement via Eq. despite 
the fact that such states are not (and are very different from) 
states. Eq. (Kll also provides a useful nonvanishing lower 
bound to the gloBal (and partial) A-particle entanglement of 
Wei states in the interval x e (1], for any even A. A compar¬ 
ison between such a bound (with D denoting the relative en¬ 
tropy), which requires only three local measurements, and the 
true value of the relative entropy of global A-particle entan¬ 
glement for these states ll5^ . whose experimental evaluation 
would conventionally require a complete state tomography, is 
presented in Eigure|^b). 


B. Genuine multiparticle entanglement 

We now show how general analytical results for geomet¬ 
ric measures of genuine multiparticle entanglement can be 
obtained as well within our approach. The results from the 
previous section, while quite versatile, cannot provide useful 
bounds for the complete hierarchy of multiparticle entangle¬ 
ment, because Ad^ states are M-separable for all M < rA/2], 
and thus in particular biseparable for any number A of qubits. 
Therefore, to investigate genuine entanglement we consider a 
different reference set of states, specifically formed by mix¬ 
tures of GHZ states, hence incarnating archetypical represen- 
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A State 

{Ci,C2,C3l 

2/=i IgI 

[9, 

|GHZ<^>) 


V3 

{cos-‘(^),0,f} 

|GHZ*''») 

{1,1,11 

3 

{0,0,01 


(- - -1 

19’ 9’ 9 / 

5 

3 

{cos->(^),0,f} 


(x, X, 2x - 1) 

2x-h |2x- 1| 

{0,0,01 

II ic<;>) 

{1,1,11 

3 

* 

icr>> 

{1,1,11 

3 

{ 1 , 0 , 0 } 


{1,1,11 

3 

{0,0,01 

|VJ/(4)) 

{1,1,11 

3 

{0,0,01 

Qf 

{1,1,11 

3 

{0,0,01 

|GHZ<^>) 


V2 



IJI _ 2 _Z_) 

1 9 V3 ’ 9 V3’ 9 V3 J 

7 

3-S^ 

|cos-‘(4),0,|| 

* 


V2x 


icf) 


3 

2 

* 

IGHZ*'”) 

{l,-l,ll 

3 

{0,0,01 

^ ewL(^) 

{x, -X, 2x - 11 

2x-h|2x- 1| 

{0,0,01 

" icf> 

{1,-1,11 

3 

* 

% icf> 

{1,-1,11 

3 

* 

|D®> 

{1,1,-11 

3 

{0,0,01 

ef 

{-1,-1,-11 

3 

{0,0,01 

r- |GHZ<’’> 


V2 

lo - —\ 

1^’ 56’ 56] 



V2x 

lo — —! 

1^’ 56’ 56/ 

icf>> 

h-i.ii 

3 

2 

* 

IGHZ"**) 

{1,1,11 

3 

{0,0,01 

00 

jx, X, 2x - 11 

2x-h |2x- 1| 

{0,0,01 

|Cf>) 

{1,1,11 

3 

* 

^ |Df> 

{1,1,11 

3 

{0,0,01 

of 

{1,1,11 

3 

{0,0,01 


TABLE III. Applications of our framework to construct ac¬ 
cessible lower bounds on global and partial (M-inseparable) 
multiparticle entanglement (which are nonzero for any 
M > \NI2'\ when 2j|Cj| > 1), for the families of N- 
qubit states listed as follows. (i) A-qubit GHZ states (43 
|GHZ<'^>) = -L (|00 • • • 00) -I-111 • • • 11» with A > 3. (ii) A-qubit W 

states (54) ^ (|00 • • • 01)- 1 -|00 • • • 10)-H ■ • •-t |10 • • • 00» 

with N > 3. (iii) A-qubit Wei states (55] |56l e^hix) = 
x|GHZ<">><GHZ<~’| -H ^ 2f=i {Pk + Pk), where A > 4, x 6 [0,1] 
and Pk is the projector onto the binary A-qubit representation of 
whereas P,- = (iv) A-qubit linear cluster states 

|C[^*) corresponding to the A-vertex linear graph • — •-• 

1141 1381 . (v) A-qubit rectangular cluster states ICj^') corresponding 

• - • - • • • - • _ 

to the vertex ladder-type graph | | | |38|[52l. (vi) 

• — • — • • • — • 

A-qubit (symmetric) Dicke states |D^^’) = 2,- n,(|0)®^“*^ 8 |1)®*^), 

which are superpositions of all states with k qubits in the excited 
state |1) and N - k qubits in the ground state |0), with the symbol 
|ni(-)lf=i denoting all the Z = distinct permutations of O’s 
and I’s; we focus on half-excited Dicke states, given by k = A/2 
for any even A (5li|52l|5ll ED- (vii) 4-qubit singlet state 1601 

= ^[|0011>-t-|1100>-(|0101)-H|0110>-H|1001>-t|1010»/2]. 

(viii) A-qubit generalised Smolin states dHlllllSD for even 
A > 4, which are instances of Al^ states with correlation triple 
j(_l)V/2 (_i)V/2 (_i)jv/2j^ hence their entanglement quantification 
is exact. The asterisk * indicates non-permutationally invariant 
states for which the optimisation of the bounds requires different 
angles for each qubit (not reported here). Notice that in the 
table we listed mostly pure states. In general, if the triple {cyj is 
optimal for a pure A-qubit state lO*'''*), then for the mixed state 
g^^\q) = -I- -^1®“, obtained by mixing with 

white noise, one still gets nonzero lower bounds to global and partial 
entanglement for all q > Ij |Cjj, as shown in Figurej^for some 
representative examples. 


tatives of full inseparability 01121 . Any such state which 
will be referred to as a GHZ-diagonal (in short, state, 
can be written as ^ = Yji,± ^ in terms of its eigen¬ 
values p*, with the eigenvectors ± crf^)|i) 

forming a basis of A-qubit GHZ states (where de¬ 

notes the binary ordered A-qubit computational basis). The 
states have been studied in recent years as testbeds for 
multiparticle entanglement detection 0E81, and specific al¬ 
gebraic measures of genuine multiparticle entanglement such 
as the A-particle concurrence Haim and negativity II3 have 
been computed for these states. Here, we calculate exactly 
the whole class of geometric measures of genuine multiparti¬ 
cle entanglement defined by Eq. (|l|, with respect to any 
contractive and jointly convex distanceD, for states of an 
arbitrary number A of qubits. 

By applying our general framework, we can prove that, for 
every valid D, the closest biseparable state to any state 
can be found within the subset of biseparable Qi^i states (see 
Supplementary Material for detailed derivations). The latter 
subset is well characterised 0, and is formed by all, and only, 
the states with eigenvalues such that = max,_± pf < 


1 /2. We can then show that the closest biseparable state to 

an arbitrary state has maximum eigenvalue equal to 1/2, 
which allows us to solve the optimisation in the definition of 
E 2 , with respect to every valid D. We have then a closed 
formula for the geometric multiparticle entanglement of any 
state ^ with maximum eigenvalue pmax, given by 




0, Pmax < 1/2; 

gD(Pmax), Otherwise, 


( 6 ) 


where go denotes a monotonically increasing function whose 
explicit form is specific to each distance D, as reported in Ta- 
ble|IV]for typical instances. 

Let us comment on some particular results. The genuine 
multiparticle trace distance of entanglement E^' is found to 
coincide with the genuine multiparticle negativity na and 
with half the genuine multiparticle concurrence ifT^ for all 
QjS! states, thus providing the latter entanglement measures 
with an insightful geometrical interpretation on this important 
set of states. Examples of Qn states include several resources 
for quantum information processing, such as the noisy GHZ 
states and Wei states introduced in the previous section. In 
particular, for noisy GHZ states (described by a pure-state 
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FIG. 3. (a) Lower bounds to the global geometric entanglement based on infidelity for noisy versions of some Ai-qubit states (defined 
in the caption of Table |ni|, as functions of the probability q of obtaining the corresponding pure states. The non-solid lines refer to bounds 
obtained by the method of Ref. (T3l for: 4-qubit linear cluster state (green dotted), 6-qubit rectangular cluster state (red dashed), 6-qubit 
half-excited Dicke state (orange dot-dashed), 4-qubit singlet state (magenta dot-dot-dashed). The solid blue line corresponds to our bound 
based on Al^-fication for all the considered states, which is accessible by measuring only the three correlation functions Cj = ((^^ 

Relative entropy of multiparticle entanglement of A^-qubit Wei states defined in Table |lll[ as a function of the probability x of obtaining 
a GHZ state. The dashed red line = x denotes the exact value of the global relative entropy of entanglement as computed in 1561 . 

The solid blue line denotes our accessible lower bound, obtained by combining Eqs. l[^ and l|^ with the expressions in Tables [n| an ding and 
given explicitly by “ log2(2 - 2x) -l- x(log 2 (x) - log 2 (l - x)) for ^ < x < 1, while it vanishes for 0 < x < The bound becomes 

tight for X = 1, thus quantifying exactly the global multiparticle entanglement of pure GHZ states. We further show that our lower bound to 
global entanglement coincides with the exact genuine multiparticle entanglement of Wei states, ^^at is computed in 

the next section of this paper. The results are scale-invariant and hold for any even N. 


Distance D 

^D(pmax) 

Relative 
entropy Dre 

1 + Pmax l0g2 Pmax 
+ (1 -Pmax)l0g2(l “Pmax) 

Trace Dxr 

pmax ~ 2 

Infidelity Dp 

2 ~ '^Pmax(l ~ pmax) 

Squared 
Bures Dp 

2 - V2 ( y 1 - Pmax + VFmax) 

Squared 
Hellinger Dh 

2 ~ V 2 ^ 1 — pmax '^pmax^ 


TABLE IV. Analytical expression of genuine multiparticle entan¬ 
glement E® for GHZ-diagonal states of any number N of qubits as 
defined by Eq. for representative choices of the distance function 
D (introduced in Table [lit. 

probability q as detailed in Table 0, we recover that every 
geometric measure of genuine multiparticle entanglement is 
nonzero if and only if q > (1 -H (1 - 2^)~*)/2 1/2 ll5l 

and monotonically increasing with q, as expected; for q — I 
(pure GHZ states), genuine and global entanglement coin¬ 
cide, i.e. the hierarchy of Figure [^collapses, meaning that all 
the entanglement of A-qubit GHZ states is genuinely shared 
among all the N particles 0^ . On the other hand, the rel¬ 
ative entropy of genuine multiparticle entanglement of Wei 
states ll55l |56l can be calculated exactly via Eq. (j^; interest¬ 
ingly, for even N it is found to coincide with the lower bound 
to their global entanglement that we had obtained by Af^- 
hcation, plotted as a solid line in Figure |^b). This means that 
for these states also the genuine multiparticle entanglement 


can be quantified entirely by measuring the three canonical 
correlation functions (c/), for any N. More generally, for arbi¬ 
trary 0!^ states, all the genuine entanglement measures given 
by Eq. ^ can be obtained by measuring the maximum GHZ 
overlap p^ax, which requires N + 1 local measurement set¬ 
tings given explicitly in Ref. Il48l . This is remarkable, since 
with the same experimental effort needed to detect full insep¬ 
arability 0 we have now a complete quantitative picture of 
genuine entanglement in these states based on any geometric 
measure, agreeing with and extending the findings of lfT4l[T8l . 
Furthermore, as evident from Eq. all the geometric mea¬ 
sures are monotonic functions of each other: our analysis thus 
reveals that there is a unique ordering of genuinely entangled 
states within the distance-based approach of Fig. 

In the same spirit as the previous section, and in compliance 
with our general framework, we note that the exact results ob¬ 
tained for the particular reference family of states provide 
quantitative lower bounds to the genuine entanglement of gen¬ 
eral A-qubit states. This follows from the observation that any 
A-qubit state g can be transformed into a state with eigen¬ 
values pf - (JSf\g\/3f} by a LOCC procedure that we may 
call GHZ-diagonalisation lfT4l . Therefore, given a completely 
general state g, one only needs to measure its overlap with a 
suitable reference GHZ state; if this overlap is found larger 
than 1/2, then by using Eq. (|^ with equal to the mea¬ 
sured overlap one obtains analytical lower bounds to the gen¬ 
uine multiparticle entanglement of g with respect to any 
desired distance D. As before, the bounds can be optimised 
in situations of partial prior knowledge, e.g. by applying lo¬ 
cal unitaries on each qubit before the GHZ-diagonalisation, 
which has the effect of maximising the overlap with a chosen 
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particular GHZ vector in the basis {(/Sf)). The bounds then re¬ 
main accessible for any state by -H 1 local measurements 
SSl . with exactly the same demand as for just witnessing en¬ 
tanglement Q. 

For instance, for the singlet state ll60l . which is 

a relevant resource in a number of quantum protocols in¬ 
cluding multiuser secret sharing one has pmax = 

^+|vj/( 4 ))^v{/( 4 )^+) _ 2/3 > 1 / 2 ^ obtainable by measuring the 

overlap with the GHZ basis state ^ 3 ) = (|0011)-i-|1100))/ V2. 
Optimised bounds to the genuine multiparticle entanglement 
of half-excited Dicke states |Dj^ 2 ) (for even N > 4), de¬ 
fined in Table ll5^ 1^ . can be found as well based on 
GHZ-diagonalisation, and are expressed by = (j^/ 2 ) 2* 
meaning that they become looser with increasing N and stay 
nonzero only up to = 8 . In this respect, we note that al¬ 
ternative methods to detect full inseparability of Dicke states 
for any N are available na i5i] EH, but quantitative results 
are lacking in general. Nevertheless, applying our general ap¬ 
proach to an alternative reference family more tailored to the 
Dicke states could yield tighter lower bounds that do not van¬ 
ish beyond N = 8 . 

Finally, notice that a lower bound to a distance-based mea¬ 
sure of genuine multiparticle entanglement, as derived in this 
section, is automatically also a lower bound to corresponding 
measures of global and any form of partial entanglement, as 
evident by looking at the geometric picture in Figure How¬ 
ever, for states which are entangled yet not genuinely entan¬ 
gled, the simple bound from the previous section remains in¬ 
strumental to assess their inseparability with minimum effort. 

states are themselves instances of such states (in fact, for 
even N, states are also Qn states, but with pmax <1/2 for 
N>2). 


C. Applications to experimental states 

In this section, we benchmark the applicability of our re¬ 
sults to real data from recent experiments Il44l 149145^ IMl . 

In Refs. IHlIia, the authors used quantum optical setups 
to prepare an instance of a bound entangled four-qubit state, 
known as Smolin state ED. Such a state cannot be written 
as a convex mixture of product states of the four qubits, yet 
no entanglement can be distilled out of it, thus incarnating the 
irreversibility in entanglement manipulation while still repre¬ 
senting a useful resource for information locking and quantum 
secret sharing ll^ l46ll . It turns out that noisy Smolin states are 
particular types of states (for any even N) Il43ll46l . that in 
the representation of Figurej^b) are located along the segment 
connecting the tetrahedron vertex {(- 1 )^^^, (- 1 )^^^, (- 1 )^^^) 
with the origin. Therefore, this work provides exact analytical 
formulae for all the nontrivial hierarchy of their global and 
partial entanglement, as mentioned in Table [nl| In the spe¬ 
cific experimental implementation of Ref 1441 for N = 4, the 
global entanglement was detected (but not quantified) via a 
witness constructed by measuring precisely the three correla¬ 
tion functions {c^). Based on the existing data alone (and with¬ 
out assuming that the produced state is within the Ad^ family). 


we can then provide a quantitative estimate to the multiparti¬ 
cle entanglement of this experimental bound entangled state 
in terms of any geometric measure E^, by using Table [h] The 
results are reported in Table |V|a) for the illustrative case of 
the trace distance. 

Remaining within the domain of quantum optics, recently 
two laboratories reported the creation of six-photon Dicke 
states |Df> lEDEa. Dicke states l58l are valuable resources 
for quantum metrology, computation, and networked com¬ 
munication, and emerge naturally in many-body systems as 
ground states of the isotropic Lipkin-Meshkov-Glick model 
ll59]l . Based on the values of the three correlation functions 
(cy), which were measured in Refs. El ED to construct some 
entanglement witnesses, we can provide quantitative bounds 
to their global and partial geometric entanglement (for 
4 < M < 6 ) from Eq. Q; see Table [V|a). 

A series of experiments at Innsbruck ll40l l49l l50l l68ll re¬ 
sulted in the generation of a variety of relevant multi-qubit 
states with trapped ion setups, for explorations of fundamental 
science and for the implementation of quantum protocols. In 
those realisations, data acquisition and processing for the pur¬ 
pose of entanglement verification was often a more demand¬ 
ing task than running the experiment itself B9ll . Focusing hrst 
on global and partial entanglement, we obtained full datasets 
for experimental density matrices corresponding to particu¬ 
larly noisy GHZ and W states of up to four qubits, produced 
during laboratory test runs Il68l . Despite the relatively low fi¬ 
delity with their ideal target states, we still obtain meaningful 
quantitative bounds from Eq. The results are compactly 
presented in Table [V|a). 

Regarding now genuine multiparticle entanglement, the au¬ 
thors of Ref. IEQI reported the creation of (noisy) GHZ states 
of up to A = 14 trapped ions. In each of these states, full 
inseparability was witnessed by measuring precisely the max¬ 
imum overlap pmax with a reference pure GHZ state, without 
the need for complete state tomography. Thanks to Eq. El- 
we can now use the same data to obtain a full quantification of 
the genuine A-particle entanglement of these realistic states, 
according to any measure E^, at no extra cost in terms of 
experimental or computational resources. The results are in 
Table|V|b), for all the representative choices of distances enu¬ 
merated in Table |lV] Notice that we do not need to assume 
that the experimentally produced states are in the set: the 
obtained results can be still safely regarded as lower bounds. 


III. DISCUSSION 

We have introduced a general framework for estimating and 
quantifying geometric entanglement monotones. This enabled 
us to achieve a compendium of exact results on the quantifi¬ 
cation of general distance-based measures of (global, partial, 
and genuine) multiparticle entanglement in some pivotal refer¬ 
ence families of A-qubit mixed states. In turn, these results al¬ 
lowed us to establish faithful lower bounds to various forms of 
multiparticle entanglement for arbitrary states, accessible by 
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a 


Global and partial multipaiticle entanglement 


State Ref. 

Fidelity (%) 


{Cl,C2,C3) 

2/=i Icjj 


gf’ ED 

96.83 ± 0.05 

{0.401 ± 0.004,0.362 ± 0.004,0.397 ± 0.008) 

1.16 ±0.01 

0.040 ± 0.002 

4 ' m 

4 Ell 

56 ±2 

{0.8 ±0.2,0.5 ±0.2, -0.3 ±0.1) 

1.6 ±0.3 

0.15 ±0.08 

65 ±2 

{0.63 ± 0.02,0.63 ± 0.02, -0.42 ± 0.02) 

1.69 ±0.04 

0.17 ±0.01 

t?(3HZ ESI 

87.9 

{-0.497,0.515,-0.341) 

1.35 

0.102 

Pghz ES] 

80.3 

{0.663,0.683,0.901) 

2.25 

0.312 

El 
4 El 

19.4 

{-0.404,0.454, -0.378) 

1.24 

0.0589 

31.4 

{0.472, 

-0.468, -0.446) 

1.39 

0.0963 

b 


Genuine multiparticle entanglement 


State Ref. 

Fidelity (%) 

^7. 

F^Tr 

^7 

^7 

^7 

Pghz El 

97.0 ± 0.3 

0.81 ±0.02 

0.470 ± 0.003 

0.329 ± 0.008 

0.36 ±0.01 

t?(3HZ El 

95.7 ±0.3 

0.74 ±0.01 

0.457 ± 0.003 

0.297 ± 0.007 

0.323 ± 0.008 

Pghz El 

94.4 ± 0.5 

0.69 ± 0.02 

0.444 ± 0.005 

0.27 ±0.01 

0.29 ± 0.01 

Pg^z El 

89.2 ± 0.4 

0.51 ±0.01 

0.392 ± 0.004 

0.190 ±0.005 

0.200 ± 0.006 

Pghz El 

81.7 ±0.4 

0.313 ±0.009 

0.317 ±0.004 

0.113 ±0.003 

0.117 ±0.003 

m 

62.6 ± 0.6 

0.046 ± 0.004 

0.126 ±0.006 

0.016 ±0.002 

0.016 ± 0.002 

El 

50.8 ± 0.9 

0.0002 ± 0.0004 

0.008 ± 0.009 

0.0001 ± 0.0001 

0.0001 ±0.0001 


TABLE V. (a) Accessible lower bounds to global and partial multiparticle entanglement of some experimentally prepared states, given by 
Eq.§ and evaluated in particular for the trace distance of entanglement by using Eq. ij^ for even N and Eq. for odd N. Following 
the theoretical analysis of Table|III| data obtained by direct measurements of the canonical correlation functions were used to construct bounds 
for a noisy Smolin state of 4 photons m, noisy Dicke states of 6 photons ElllsU, and noisy GHZ states of 4 ions (68]. For noisy GHZ 
states of 3 ions and noisy W states of 4 ions (two implementations labelled as A and B) l68l . full datasets were used to extract the optimised 
correlation functions (cy) required for the bounds. For all the presented experimental states (whose fidelities with the ideal target states are 
reported for reference), we are able to provide a reliable estimate of the multiparticle entanglement for any M > \N 12]. (b) Lower bounds 
to genuine multiparticle entanglement of expe rime ntal noisy GHZ states of up to 14 ions (50), as quantified in terms of all the distance-based 
entanglement measures E® reported in Table IV obtained by Eq. with given in each case by the measured fidelity with the pure 
reference GHZ state. All the reported entanglement estimates are obtained from the same data needed to witness full inseparability, which for 
general A-qubit states can be accessed by A -l- 1 local measurements without the need for a full tomography. 


few local measurements and effective on prominent resource 
states for quantum information processing. 

Our results can be regarded as realising simple yet par¬ 
ticularly convenient instances of quantitative entanglement 
witnesses 1221123], with the crucial advance that our lower 
bounds are analytical (in contrast to conventional numerical 
approaches requiring semidefinite programming) and hold for 
all valid geometric measures of entanglement, which are en¬ 
dowed with meaningful operational interpretations yet have 
been traditionally hard to evaluate lfT3ll69ll . 

A key aspect of our analysis lies in fact in the generality 
of the adopted techniques, which rely on natural information- 
theoretic requirements of contractivity and joint convexity of 
any valid distance D entering Eq. Q- We can expect our gen¬ 
eral framework to be applicable to other reference families of 
states (for example, states diagonal in a basis of cluster states 
011691, or more general states with X-shaped density matri¬ 
ces m), thereby leading to alternative entanglement bounds 
for arbitrary states, which might be more tailored to different 
classes, or to specific measurement settings in laboratory. 

Furthermore, our framework lends itself to numerous other 
applications. These include the obtention of accessible ana¬ 
lytical results for the geometric quantification of other useful 
forms of multiparticle quantum corTelations, such as Einstein- 
Podolsky-Rosen steering Gomii, and Bell nonlocality in 
many-body systems 15^ . This can eventually lead to a uni¬ 


fying characterisation, resting on the structure of information 
geometry, of the whole spectrum of genuine signatures of 
quantumness in cooperative phenomena. We plan to extend 
our approach in this sense in subsequent works. 


Another key feature of our results is the experimental acces¬ 
sibility. Having tested our entanglement bounds on a selection 
of very different families of theoretical and experimentally 
produced states with high levels of noise, we can certify their 
usefulness in realistic scenarios. We recall that, for instance, 
three canonical local measurements suffice to quantify exactly 
the global entanglement of GHZ states of any even number N 
of qubits, while N + 1 local measurements provide their exact 
genuine entanglement, according to every geometric measure 
for any N, when such states are realistically mixed with white 
noise. Compared to other complementary studies of accessi¬ 
ble quantification of multiparticle entanglement lfT3lj23] . our 
study retains not only a comparably low resource demand but 
also crucial aspects such as efficiency and versatility, as shown 
in Table [I] This can lead to a considerable simplification of 
quantitative resource assessment in future experiments based 
on large-scale entangled registers, involving e.g. two quantum 
bytes (16 qubits) and beyond Il50ll68ll . 
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IV. METHODS 


Distance-based measures of multiparticle entanglement. 

A general distance-based measure of multiparticle entangle¬ 
ment is defined in Eq. ([l|. In this work, the distance 
D is required to satisfy the following two physical constraints: 
(D.i) Contractivity under quantum channels, i.e. 
D(Q.{q),Q.(q')) < D{g,g'), for any states g, g', 

and any completely positive trace preserving map Q; 
(D.ii) Joint convexity, i.e. D(qg H- (1 - q)g', qx + i} — q)x') ^ 
qD{g,x) + (1 - q)D{g',x')^ for states g, g',x^ ^^d 
x', and any q e [ 0 , 1 ]. 

Constraint (D.i) implies that is invariant under local uni- 
taries and monotonically nonincreasing under LOCC (i.e., it 
is an entanglement monotone ||24l). Constraint (D.ii) implies 
that is also convex. A selection of distance functionals 
respecting these properties is given in Table [II] 


Al^-flcation. Theorem. Any A-qubit state g can be trans¬ 
formed into a corresponding AI^ state m through a fixed 
transformation, 0, consisting of single-qubit LOCC, such that 
©(£>) = ttr = ^ (l®'^ + ZLi c;cr®^), where a = Tr(po-®'^). 
Proof. Here we sketch the form of the Al^-fication chan¬ 
nel 0. We begin by setting 2{N - 1) single-qubit local 
unitaries = {(cri ® cri igi I®^'^), (I ® cri ® cri ® 

F^-3), . . . , (I®^-3 ® 0-1 ® 0-1 ® I), (I®'^-^ ® o-j ® cTi), ( 0-2 ^ 
(I ® 0-2 ® crz ® I®^-3), . . . , 


0-2 


' cr2 Qs) cr2 ' 


I), (I®^ ^ ® 0-2 ® 0 - 2 )}. Then, we fix a sequence of states 
{po,Pi,.. .p 2 (iv-i)) defined by gj = f (p^.i -H t/yp^_it/j) , for 
j e {1,2,. ..2(N - 1)). By setting go = g and p 2 (v-i) = ©(p), 

'^'""U'guf, 

where: 


we define the required channel: 0 (^) = ^ 2 (a-i) 

^22(iV-l) _ 

=1 ” 




[U, 


2(N-l) 


l'l2>ll = l’ ■ ■ ■ ’ 

,Notice that is still 

a sequence of single-qubit local unitaries. Since 0 is a con¬ 
vex mixture of such local unitaries, it belongs to the class of 
single-qubit LOCC, mapping any M-separable set into itself. 
In the Supplementary Material, we show that 0(p) = vi, con¬ 
cluding the proof 


Lower bound optimisation. Lor any valid distance-based 
measure of global and partial multiparticle entanglement E^, 
the maximisation in Eq. (j^ is equivalent (for even N) to max¬ 
imising |ci| + \c 2 \ + IC 3 I, where cj - Tr[t/gipI/gCr®^], over 
local single-qubit unitaries (a = 1 ,...,A). 

By using the well known correspondence between the spe¬ 
cial unitary group SU(2) and special orthogonal group S0(3), 
we have that to any one-qubit unitary corresponds the or¬ 
thogonal 3x3 matrix 0^“' such that = ((9*“^n)-(^, 

where ft = { 111 , 112 , 112 } e and cf = (cri, 0 - 2 , 0 - 3 ) is the 
vector of Pauli matrices. We have then that sup|(j(„)|(|ci| -H 
IqI -I- IC 3 I) = sup|Q(„)|(|fii...i| + |f 22 - 2 l -I- IL 33 ... 3 I), where 
T. . • — V T ■ ■ ■ nnH t. . . — 


Tr [p (cr,j ® cT/j ® ■ ■ ■ ® cr, 2 ^)]. In the case of permutationally 
invariant states g, the 3 x 3 x ■ ■ ■ x 3 tensor is fully 

symmetric, i.e. for any permutation 1 ? 

of the indices, so that the optimisation can be achieved when 
(9^1 = - ■ ■ ■ - fjW iIt^ . As indicated in the main text, 

we then need to perform the maximisation over just the three 
angles ( 0 , tjj, f] which determine the orthogonal matrix 0 *“^ 
corresponding to an arbitrary single-qubit unitary 




Q 

COS ' 2 

. . Q : ±:± 

-I sin je 2 


e ^ 


-i- 


-i sin 

f) i'tii. 

cos 2 


\ 


/ 


As a special case, for a two-qubit state (N - 2) the optimal 
local operation is the one which diagonalises the correlation 
matrix 
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SUPPLEMENTARY MATERIAL 
Appendix A: Multiparticle entanglement 

When considering a multiparticle quantum system, there 
exist two different approaches to entanglement, one referring 
to a particular partition of the composite system under consid¬ 
eration (partition-dependent setting), and another which con¬ 
siders indiscriminately all the partitions with a set number of 
patties (partition-independent setting). 

In order to characterise the possible partitions of an A-qubit 
system, we will employ the following notation ll32l : 

• the positive integer M, I < M < N, representing the 
number of subsystems; 

• the sequence of positive integers := 

{Ki,K 2 ,"-,Km}, where a given Ka represents 
the number of qubits belonging to the cr-th subsystem; 

• the sequence of sequences of positive integers 

such that Qa - {/j“\ • • • , with e {1, • • ■ ,N} 

and Qa (T Qa’ - 0 for a + a', where a given sequence 
Qa represents precisely the qubits belonging to the a-th 
subsystem. 

In the following we will say that {2 q.)^j identifies a generic 
M-partition of an A-qubit system. 

The set of A-qubit separable states ^ with respect to 

the M-partition {Qa]a=\ contains all, and only, states g of the 
form 

^ ® T® ® . . . ® (Al) 

i 

where {/?,) forms a probability distribution and are ar¬ 
bitrary states of the a-th subsystem. In other words, any 
{2Q,)^j-separable state can be written as a convex combina¬ 
tion of product states that are all factorised with respect to the 
same partition [Qd^j. On the other hand, the set of A-qubit 
M-separable states Sm contains all, and only, states that can 
be written as convex combinations of product states, each of 
which is factorised with respect to an M-partition that need not 
be the same. One can easily see that the set of M-separable 
states is the convex hull of the union of all the sets of {2a)^j- 
separable states obtained by considering all the possible M- 
paititions (Qa)" j. 

Any valid measure of multiparticle entanglement must be 
zero on the relevant set of separable states and monotonically 
non-increasing under LOCC. In the partition-dependent set¬ 
ting, a LOCC with respect to a particular -partition 

amounts to allowing each of the M parties to perform local 
operations on their qubits, and communicate with any other 
patty via a classical channel |[3l. Conversely, in the partition- 
independent setting, one considers operations that are LOCC 
with respect to all of the M-partitions, which can be shown to 
be all and only the single-particle LOCC. A convex combina¬ 
tion of single-particle local unitaries acting on a state g, given 
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by 


and 




PiUf''®ufh 


(A2) 


d* = ^ [i + Cl + {-if^\-iyc 2 + (-1FC3], 


(B3) 


is a particular type of single-particle LOCC (requiring only 
one-way communication). It can be physically achieved by 
allowing one of the subsystems a to randomly select a local 
unitary t/f“^ by using the probability distribution (p,) and then 
to communicate the result to all the other subsystems. 

One can also impose that a measure of A-particle entangle¬ 
ment is convex under convex combinations of quantum states, 
i.e. 


“ ^)e') ^ “ ?)^IG„)", 

Eniqg + (1 - q)Q') < qEuig) + (l - q)EM(g') (A3) 


for some probability q and quantum states g and g', which en¬ 
sures that classical mixing of quantum states cannot lead to 
an increasing of entanglement. Any measure obeying these 
properties is referred to as a convex entanglement monotone. 
In this work we adopt geometric measures of multiparticle en¬ 
tanglement, defined in terms of the distance to the relevant set 
of separable states. For a given distance D, generic distance- 
based measures of the multiparticle entanglement of an N- 
qubit state g, quantifying how much g is not { 2 q.)“ j-separable 
(resp., M-separable), are given by, respectively, 

« “-1 

E^(g)= inf D(g,q), (A5) 


where Eq. (A4i refers to the partition-dependent setting, and 
Eq. ( |A5| l to the partition-independent one. It is sufficient for 
the distance D to obey contractivity and joint convexity (see 
Methods in the main text) for £; 
entanglement monotones l24l . 


IQaC 


and to be convex 


where i 6 { 1 , ■ ■ ■ , 2 ^ '), is the binary ordered A-qubit 

computational basis and finally p is the parity of j/Sf) with 
respect to the parity operator along the z-axis 113 = cr|^, i.e. 

n3l6f) = (-iyi/3f). (B4) 


In the {ci,C 2 ,C 3 )-space, the set of states with even 
N is represented by the tetrahedron 7 '(_i)iv /2 with ver¬ 
tices { 1 ,(- 1 )^/^ 1 ), {-l,-(-l)'^'' 2 , 1 ), {l,-(-l)^/^-l) and 
as illustrated in Eig. 2(b) in the main text. 
This tetrahedron is constructed si mply by imposing the non¬ 
negativity of the four eigenvalues (B3 1 of such states. 


Eor odd N, the eigenvectors and eigenvalues of the AI^ 
states can be easily written in spherical coordinates as 


laf) - cos 


- + (i + (-iy)- 
2 4J 


10 


(B5) 


+(-iye' 




6 „ TT 


mN\ 


10 , 


and 


4, = ^(l±r), (B6) 

where i 6 (I,-- - ,2^“*), {lOl^^i again the binary ordered 
A-qubit computational basis, p is the parity of |i) with re¬ 
spect to the parity operator 113 = cr|^, ci = rsin 0 cos(^, 

C 2 = r sin 6 sin 0 and — r cos 0, with r = + <22 + 

6 6 [ 0 , tt] and (p 6 [0,27r[. 

Consequently, thanks again to the semi-positivity con¬ 
straint, the set of A1^ states with odd A is represented in the 
{ci, C 2 , C 3 )-space by the unit ball 'B\ centred into the origin, as 
shown in Eig. 2(c) in the main text. 


Appendix B: The set of Af^ states 


In this appendix we show some relevant properties of the 
subclass of A-qubit states with all maximally mixed marginals 
that we refer to as AI^ states. Their matrix representation in 
the computational basis is the following; 


m - 


1 



CiCr 


.»N 


(Bl) 


where I is the 2x2 identity matrix, cr,- is the i-th Pauli matrix, 
Ci - Tr |^uTcr®^j G [-1,1] and A > 1. These states are denoted 
by the triple {ci,C 2 ,C 3 ). 

The characterisation of the AI^ states is manifestly differ¬ 
ent between the even and odd A case. Eor even A, the eigen¬ 
vectors and eigenvalues are given by, respectively. 


Appendix C: A4^-fication 


The following Theorem is crucial for providing a lower 
bound to any multiparticle entanglement monotone of any 
state g and for analytically computing the multiparticle geo¬ 
metric entanglement of any Mjf state m. 

Theorem C.l. Any N-qubit state g can be transformed into 
a corresponding A4^ state through a fixed operation, 0 , 
that is a single-qubit LOCC and such that 


0(e) = gmi = 


1 

2A 



(Cl) 


where ct - Tr(pcr®^). 


\Pt) = 


1 



±Oio, 


(B2) 


Proof The first part of the proof was sketched in the Methods 
section and is repeated here for completeness. 
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We will give the form of ©(£)), show that 0 is a single-qubit 
LOCC, and finally prove that it transforms any A^-qubit state 

pintopyvt^. 

To define ©(p), we begin by setting 2{N - 1) single-qubit 
local unitaries 

- {(‘^1 ® o"! ® 7®^“^), (7 (g) cTi ® cTi ® 7®^“^), 

... (7®^~3 ® O’! ® cTi ® 7), (7®^“2 ® cTi ® (Ti) 

, (cr2 ® 0-2 ® (7 ® 0-2 ® 0-2 ® 7®^“^), 

. . . (7®^~3 ® 0-2 ® 0-2 ® I), (7®^“^ ® 0-2 ® 0 - 2 )). 

(C2) 


is clear that the R^'. . of oi aie . . if q and h are (i) any 
combination of only 1 and 0 or (ii) any combination of only 2 
and 3, and zero otherwise. 

Generally, for j € [1,77 - 1], the 7?^'. . of p, are R‘:^t' . 
if ij and ij+\ are (i) any combination of only 1 and 0 or (ii) 
any combination of only 2 and 3, and zero otherwise. For j e 
[77,2(77 - 1)] the conditions are analogous, where the 7?^7^ 
of Qj are if ij and ij+i are (i) any combination of only 

2 and 0 or (ii) any combination of only 1 and 3, and zero 
otherwise. For the final state p 2 (v-i). the only nonzero 
are those for which (/ 1/2 ... in} consists of only 0, 1,2, or 3, 
and that for these elements - R^^. . . Therefore 


Then, we fix a sequence of states (po,Pi, ■ ■ ■ 7 ' 2 (v-i)) defined 
by 


Qj^\{ej-i+Ujgj-iU]) (C3) 

for ;■ e (1,2,... 2(77 - 1)). By setting go = g and p 2 (v-i) = 
©(p), we define the required LOCC channel, i.e. ©(p) = 
22 (tiy ToJi ' U[gU'.^ where U'. are the following unitaries 


J 3 

©(e) = e 2 (iv-i) = ^ ^ 


2N 

1 

7^ 


:(Ti (SCTi . . .is CTi 


!=0 

/ 




i=l 


- eAii (C8) 


where we have used R^. . - Tr(pcr®^) s c, for / e (1,2,3) and 
■^ 00...0 “ Tr(p) = 1 . 


{u']tr - 


{Ui 

{Ui.Uo) 


,2(V-1) 
hi = l 
2(V-1) 
i2>i\=l 


i 2 (W-l) 


{^I2(JV-1) ■ ■ ■ ^'2 ^<1 );2(«-l)>...>i2>li=l 


(C4) 


It is clear that are unitaries that still act locally on 

individual qubits. Since © is a convex mixture of such local 
unitaries, we conclude that © is a single-qubit LOCC. 

Now we will show that ©(p) = gM},- Consider the arbitrary 
77-qubit state p written in the form 


ON Zj 


2N /_i 


)cri. 


(C5) 


where the = Tr [p cr,j ® cri^... ® cr,-^] e [-1,1] are 

the correlation tensor elements of p with ctq = I. Convex 
combination of two arbitrary A^-qubit states p and p' gives 

, 3 

pp -H (1 - q)g' ^ ^ Yj ®’ c ^'-2 ■ ■ ■ ® o-iN (C 6 ) 

/ 1 / 2 ... 177=0 


Herein, we will refer to pyvt^ = ©(p) as the Af^-fication of 
the state p. Theorem |C.l| has two major implications. The first 
implication applies to any multiparticle entanglement mono¬ 
tone, be it partition-dependent or independent. We have that 

EmigMi) ~ ^M(®ig)) < Emig), (CIO) 

where in the first equality we use pyvr^ = ©(p) and in the in¬ 
equality we use the monotonicity under single-qubit LOCC of 
any measure of multiparticle entanglement and the fact that © 
is a single-qubit LOCC. In other words, the multiparticle en¬ 
tanglement of the Al^-fication p^y,^ of any state p provides us 
with a lower bound of the multiparticle entanglement of p. 

The second implication applies specifically to distance- 
based measures of multiparticle entanglement, although re¬ 
gardless of whether such a measure is partition-dependent or 
independent. We have that, for any state m and any sep¬ 
arable state g, 

D(zu, g^i) - D(&(zu), 0(p)) < D(m, g), (Cl 1) 


where - qR^. . -i-(l-p)7?^. . . 

We will now understand the evolution of the R^:^. . for 

_ __ 

each step j in Eq. (|C3|l. The action of Ui on p is 




2 W 2 J 

M^2---^77=0 

®cr,- 


) cricr ucr 1 


(C7) 


From cricr,cri = -(-l)'^"'’^^‘'cr,- we have that the cor¬ 


relation tensor elements of U\gU[ are R^ 

(- 1 )^ 0 ,,+c5i,[+ 50 , 2 + 51 . 277 ? . . By using Eq. (C3 1 and Eq. (|C 6 |), it 


iii2...iN 


where in the first equality we use the invariance of any AI^ 
state through © and that ©(p) s is the Al^-fication of 
g, and in the inequality we use the contractivity of the dis¬ 
tance through any completely positive trace-preserving chan¬ 
nel. Moreover, the Ad^-fication of any separable state g, 
regardless of whether g is -separable or M-separable, 

is a separable Ad^ state of the same kind as g, since © is a 
single-qubit LOCC and thus leaves any set of separable states 

invariant. Therefore, both the sets >S A „ and of, respec- 

\\2a]Q-\ 

tively, ((2a)^i-separable and M-separable Ad^ states will be 
crucial to identify (see Appendix [D]i, since they allow us to 
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use Eq. ( Cl l| l to say that for any distance-based measure of 
multiparticle entanglement of an state m. 


"(Ga 


(m) = inf 

1 


D(m, g) = 


inf 

M 




E^{m) = inf Dim, g) - inf 


?<^Sm 


D(m,gMl\ 

ISul" 1 

D{Tn,gf^), 


(C12) 


^ m3 


i.e. that one of the closest {Q^l^i-separable (resp., M- 
separable) states g^ to an state ur is itself an Af^ state. 
We now formalise these two results as corollaries. 


Corollary C.l. For any N-qubit state g, the multiparticle en¬ 
tanglement of the corresponding M^j^-fied state is always 
less than or equal to the multiparticle entanglement of g, i.e. 

'^(G„)"/£'AtJ,) ^(C13) 
Em{Qmi) - (C14) 

for any {Qfi’^^^-partition of the N-qubit system and any 2 < 
M <N. 

Corollary C.l. For any contractive distance D and any Ad^ 
state vj, one of the closest {Q„]^^^-separable (resp., M- 
separable) states gj^ to m is itself an Ad^ state, i.e. 




1 

fN 




JSN 


(C15) 


for any {Qa]’^^^-partition of the N-qubit system and any 2 < 
M<N. 

Theorem IC. H allows for another result which will be useful 
to characterise the set of separable Ad^ states. 

Corollary C.3. The set of the triples {ci,C 2 ,C 3 ), with Ci — 
Tr(pcr®^), obtained by considering any possible N-qubit state 
g is 


• the unit ball Si, when N is odd; 

• the tetrahedron when N is even. 

This is because the set of Ad^-fications of all the states co¬ 
incides exactly with the set of Ad^ states. Indeed, the Ad^- 
fication channel 0 makes the entire set of states collapse into 
the set of Ad^ states, whereas it leaves the set of Ad^ states 
invariant. Herein, we shall refer to the triple {ci,C 2 ,C 3 ), with 
Ci - Tr(pcr®^), as the Pauli correlation vector corresponding 
to the state g. 


Appendix D: The set of separable Ad^ states 

We are now ready to characterise the sets and 

of, respectively, {2Q,)^j-separable and M-separable Ad^ 
states. The first ingredient is to note that S " ^ coincides 

ItiO'/Q,—1 


exactly with the set j of the Ad^-fications of any 

{fiffl^i-separable state. Furthermore, we note that since any 
Ad^ state is invariant under any permutation of the N qubits, 

then the set of {2cr)“,^i-separable Ad^ states does not 

depend on which qubits belong to each of the subsystems. 
Therefore we need only to specify the cardinalities {Ka]^^y 

to completely characterise S ^ ^ , and we will herein refer to 
the latter as the set of .-separable AdT states >5 % ■ 

Theorem D.l. For any N, the set of separable Ad^ states 
Af 

•S, f.M is either 

• the set of all Ad^ states, for any allowed {Ka]^^i parti¬ 
tion such that Ka is odd for at most one value of a; 

• the set of Ai^/ states represented in the {ci,C 2 ,C 3 )- 

space by the unit octahedron Oi with vertices {+1,0,0), 
{0,±1,0) and {0,0,±1), for any allowed par¬ 

tition such that Ka is odd for more than one value of 

a. 


Proof 

In order to characterise the set of -separable Ad^ 

states, <S, „ , we simply need to identify its representation in 

Ai^ r 1 

the {ci,C 2 ,C 3 {-space. Since = 0 J, we know 

that such a representation is the set of Pauli correlation vectors 
corresponding to all the elements of ^. 

Due to Eq. the Pauli correlation vector of any g € 

given by 


Sj = Tr = Tr ^ ® t\ 


JM) 




Ti CTj 


>T. CTj 


>T cr. 


M 


M 


=Z/’'n’'3W"+r-)=i;p'n'S’ 


(Dl) 


i Q'=l 


/ (T=l 


where in the final equality we denote c^‘^j - Tr^T|“V®^“) 
as the j-th component of the Pauli correlation vector = 
{C;“\c| 2 \ corresponding to the arbitrary state of sub¬ 
system a. Eq. (Dll can be simplified further by introducing 


the Hadamard product as the componentwise multiplication 
of vectors, i.e. for m = {ui,U 2 ,u^} and v = {vi,V 2 ,V 3 ) the 
Hadamard product is m o it = {miVi,M2V2,M3V3). Using the 


Hadamard product gives Eq. (Dl 1 as 




Al) A 2 ) 
cr. o cr. ' 


-4M) 

o cr. r 


(D2) 


i.e., that the Pauli correlation vector of any {Qal^j-separable 
state is a convex combination of Hadamard products of Pauli 
correlatio n vec tors corresponding to subsystem states. Due to 
Corollary C.3 we know that cf“^ 6 Si when Ka is odd and 
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€ T( 


Mi, 




when Ka is even, and so is represented 


by the following set 


{K„Y 


conv 


with 


^(„) ^ I Si if Ka is odd, 
T(^\)Kaii if Ka is even. 


(D3) 


(D4) 


where we define the Hadamard product between any two sets 
A and B as A o B = {a o b \ a e A, b e B} and the con¬ 
vex hull conv(A) is the set of all possible convex combina¬ 
tions of elements in A. The commutativity and associativity 
of the Hadamard product allow us to rearrange the ordering in 
Eq. ( |D3| l in the following way 

rr _ \ 


JW 

{K,. 


0^(-i)V^ ° Os 




(D5) 


where O" ° o ... o 

^0'=1 

By writing any vector in T±i as a convex combination of 
the vertices of T±i, one can easily show that 


T--ior_i =7-1, 
TioT.^Tu 
7-1 o Til = Til, 


so that 


0^(- 


(_l/p/2 - 7“(_1)A1- 


(D6) 


(D7) 


where Af- is the number of with odd Similarly, one 
can see that 


71,1 o Si = Si. 
Finally, we have that 

conv (OlliSi) - 0\ Vn > 2. 


(D8) 


(D9) 


Indeed, since {{+1,0,0), {0,+1,0), {0,0,+ 1)) c 0",Si and 
conv{{ + l, 0,0), (0, +1,0), (0,0, +1)) = (9i, we know that 0\ c 
conv(OfciSi). Now we will show that 0\ 2 conv(0”iSi). 
To do so, it is sufficient to see that 


bob' eOi 


(DIO) 


for any b,b' e Si, which trivially implies that 0"|Si £ Oi, 
andsoconv(0"iSi) £ conv{Oi) -Oi. Equation (DIO 1 holds 
since 

\bib'^\ + 1^2^21 + = l^ll l^'il + 1^2! 1^2! + l^3l 1^3! 

- n ■ ft' - ||n|| ||n'|| cos 6 < I, 

(Dll) 

where we define n = {|^iU^2U^3l) and ft' - {|l7j|, 
respectively, as the vectors corresponding to b and b' in the 
positive octant of the unit ball, and 6 as the angle between 
these vectors. 

Now, due to Eqs. (D5 1, (D71, (D8 1 and ( |D9| l, and the fact 
that conv(A) = A for any convex set A, we identify four cases; 


1. \f Ka is even for any a then 


{K„. 


Or, 




= conv 

= 71 


(T'(-I)AI-) 


(D12) 


where Ai- is the number of Kfj with odd K^jj!-, 
2. if Ka is odd for just one value of a then 


_ 


= conv 


O r,_if,,i2 1 o Si 


\p-.K^even 

— conv (T±i o Si) 
= Si; 


3. if Ka is odd for all values of a then 


(D13) 


c*'' w 

= Oi; 


Osi 


(D14) 


4. otherwise, 

_ 


O r,,,KAo\QfB, 


\fi’K^e^ 


= conv 


77,1 o (O Si 


= conv [T±\ o Si o ... o ^ J 


= conv 1 z>\ 

v-.Kyodd 


Osi 




(D15) 


For any even N-qubit system, only a {Ka}^ partitioning 
within cases[2|^and|^may be realised. In case[T] i.e. when Ka 


oAli 


is even for any a, we have S, , where M- is the 

number of Ka with odd Kal2. However, one can simply see 


that(-l)^- = (-1)^^^, and thus Sj 




M\, states. Otherwise, incases 


j^andj^ 


is the setT”(^i)W2 of all 
wehaveSj^"„^ =<9i. 

\M 


For any odd N-qubit system, only a partitioning 

within cases @[3] and may be realised. In case[^ i.e. when 


Ka is odd for only one a, we have S 




- Si, and thus 


-,Mi, 


is the set Si of all AI„ states. Otherwise, in cases 


and 


we have S 


{KaY 


= Oi. 
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By identifying the set of separable states S ^ „ , The- 
orem |D.l] implies the following corollary. 

Corollary D.l. For any multiparticle entanglement mono¬ 
tone ^ and any state mpartitioned along any given 

{Qa}^^i-partition, ^Oif 

1. Ka is odd for at most one value of a; 

2. Kg, is odd for more than one value of a and |ci| + \c 2 \ + 
kal < Iforci = Tr(zrrcr®^). 

Now we are ready to characterise also the set of M- 
-2 

separable Af^ states Indeed we know that <5^" is 

just the convex hull of the union of all the sets of {Qa]^^i- 
separable Ad^ states obtained by considering all the 

possible M-partitions {Qa]^^y Furthermore, one can easily 
see that for any M < \NI2'\ one can always find an M-partition 
such that Ka is odd for at most one value of a and thus 

'1 

whereas for any M > \NI2'\ this is impossible 
and thus „ = Oi for any possible M-partition 

I \la i 

This immediately implies the following two Corollaries. 

Corollary D.l. For any N, the set of M-separable Ad^ states 
is either 

• the set of all Ad^ states, for any M < \N 12]; 

• the set of states represented in the {ci,C2,C3)- 
space by the unit octahedron Oi with vertices {±1,0,0), 
|0, +1,0) and (0,0, ±l),/or any M > \NI2\ 

Corollary D.3. For any multiparticle entanglement monotone 
Em and any Ad^ state m, EMirn) — 0 if 

1. M < \NI2\ 

2. M> \N|2^ and |ci| + |c2| + ksl < 1/or a = Tr(nTcrf^). 


Appendix E: Multiparticle entanglement of Ad^ states 

We now provide the analytical expressions for both 
the partition-dependent and partition-independent geomet¬ 
ric measures of multiparticle entanglement E® „ (ur) and 

E^i-ur) of any Ad^ state ur. Within the partition-dependent 
setting we will restrict to any nontrivial partition i.e. 

such that Ka is odd for at least two values of a, whereas within 
the partition-independent setting we will restrict to any non 
trivial number of parties M', i.e. such that M' > [N/2']. Ac¬ 
cording to Appendix [C| and [P] in both cases we simply need 
to find the minimal distance from m to the set of Ad^ states 
inside the unit octahedron Oi. In the even N case, the closest 
state is the same for any convex and contractive distance (note 
that every jointly convex distance is convex), while in the odd 
N case this is not true. 


1. Even N case 


For even N, both the and M'-inseparable Ad^ 

states belong to the four corners obtained by removing the 
unit octahedron Oi from the whole tetrahedron 7'(_i)«/2 of Ad^ 
states. In the following we will focus only on the corner con¬ 
taining the vertex (-1, (-1)^'^^,-1), since all the Ad^ states 
belonging to the other three corners can be obtained from this 
by simply applying a single-qubit local unitary cr, ® 
i 6 (1,2,3), under which any sort of multiparticle entangle¬ 
ment is invariant. 

In order to characterise all the Ad^ states with even N be¬ 
longing to the (-1, (-1)^^^,-l)-corner, it will be convenient 
to move from the coordinate system {ci,C2,C3) to a new co¬ 
ordinate system {p,q,h), where we assign the coordinates 
Q, j, l) to the vertex {-1,(-1)^^^, -1) and the coordinates 


_ 1 + Cl - (-1)^^^C2 - C3 
^ 2 -V C\ - {-1^1^02 +cf 

_ 1 + Cl + (-1)^^^C2 + C3 
^ 2 -V Cl-{-1^1^02 +cf 

/,^(_1_(Ci-(-1)^/2c2+C3))/2, 


(El) 

(E2) 

(E3) 


to any other point in the corner. In order to avoid confusion 
between the above two coordinate systems, we will denote an 
Ad^ state m with curly brackets when representing it in the 
{ci,C 2,C3) coordinate system, whereas we will denote m with 
round brackets when representing it in the {p, q, h) coordinate 
system. Specifically, the Ad^ states represented by the triples 
{p,q,h), with a hxed value of h 6 [0,1[, correspond in the 
{ci,C2,C3)-space to all, and only, the Ad^ states belonging to 
the triangle with the following vertices: 

Vi(h) = [-h,(-lf/^h,-l], 

V2(h)^[-h,{-\fl^,-h], (E4) 

V^ih) = [-\,{-\f^^h,-h], 

in such a way that 


ip,q,h) = pVi(h) -V qV 2 (h) + {I - p - qW^ih). (E5) 


These triangles corresponding to constant values of h will play 
a crucial role, as they represent the sets of Ad^ states with con¬ 
stant {Ka]^^y- and M'-inseparable multiparticle entanglement 
for even N. In particular, for h = 0 we get one of the faces of 
the octahedron of {Ai^)^j- and M'-separable states, whereas 
with increasing h, we will prove that both the {A’q,)^j- and 
M'-inseparable multiparticle entanglement of the Ad^ states 
belonging to the corresponding triangle will increase mono- 
tonically. We will now show that the {AiQ,)“ j-separable (resp., 
M'-separable) state represented by the triple (p, q, 0) is one of 
the closest {/r^)^j-separable (resp., M'-separable) states to 
the Ad^ state (p, q, h). 

Lemma E.l. For every even N, according to any convex 
and contractive distance, one of the closest {Kfi’^^-y-separable 
(resp., M'-separable) states to any Ad^ state m belonging 
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to the {— 1 , (——l)-comer is always an state of the 
form (p',q',0) for some p',q' e [ 0 , 1 ], p' + q' < 1 . 


Proof. Let m and q be, respectively, any state belonging 
to the {- 1 , (- 1 )^''^, -l)-corner and any -separable Ad^ 

state, i.e. any Ad^ state contained in the unit octahedron 0\. 
There will always be a {-separable Ad^ state q', belong¬ 
ing to the octahedron face whose vertices are Vi ( 0 ), 1 ^ 2 ( 0 ), and 
V 3 ( 0 ) given in Eqs. (E4i, such that q' - Am+(l-A.)q for some 
d e [0,1]. Now, for any convex distance, the following holds 


D(m,q') (E 6 ) 

= D(m, Am H- (1 - /!)§•) 

< AD(m, m) + ( I — A)D(m, q) 

= (1 - A)D(m,q) 

< D(m, q). 

As one of the closest {/f^j^j-separable states q.^^ to any Ad^ 
state m is always a {/f^j^j-separable Ad^ state, then the 
above inequality implies that, for any Ad^ state m belonging 
to the {-1, (-1)^^^,-l)-corner, q^ always belongs to the tri¬ 
angle with vertices yi( 0 ), 1 ^ 2 ( 0 ), and V’ 3 ( 0 ) i.e. q^ = (p', q', 0 ) 
for some p',q' e {0,1}, p' + q' < 1. Exactly the same 
proof holds when substituting -separability with M'- 

separability. 


Lemma E.2. For every even N, any contractive distance sat¬ 
isfies the following translational invariance property: 

D (ip, /z), (p, 0)) = DII ii, /zj , J ^^, 0j j, (E7) 


which in turn is due to Eqs. (E4 1 , (E51, and: 


t/l{ci,C2,C3)t/J = {-(-1)'^/2^2,-(-1)'^^"c3,ci), (E13) 
U2{CuC2,C3}Ul ^ {C3,-(-lf/W-(-lf/^C2}. (E14) 


In order to prove the opposite inequality and thus Eq. (E7 1 , 
we now introduce a global A-qubit channel Q. with operator- 
sum representation 


2N 

^{Q)^^AiQA\, (E15) 

(=1 

where 

■!!', {I'i'jx'pjiiJl (E16) 

with the 2^ Kraus operators satisfying 2, A^A, = I®^, where 
(lO*)) and {I'Pp) constitute the binary ordered A-qubit eigen¬ 
vectors {\/3f}} with even and odd parity, respectively, i.e. they 
are such that 


= \^p, 

n3|T'7) = -|'I'7>, (E17) 


where y 6 {1, • • • ,2^ ^). It will be crucial in the following 
to see that the effect of Q on an Ad^, state represented by the 
triple Q, I, /z) is given by 


for any p,q & [ 0 , 1 ] with p -i- q < I and /z 6 [ 0,1 [. 

Proof. Eirst of all, by considering the following single-qubit 
LOCC, 



(1,0,/z). 


(E18) 


?)(£’) ^ PQ + qU\QU\ + (1 - E - q)U 2 gUl (E8) 


where p,q € [ 0 , 1 ], p H- ^ < 1 and 

Ui = (E9) 

U2 = (ElO) 

with 5,- - (I + * 0 ",) and we have 

the following inequality. 


D 




" (P^q’h),A^i_y_^{p,q,Q)) 

< D({p, q, h) , (p, q, 0)), 


(Ell) 


where the final inequality is due to the contractivity of the 
distance D, whereas the first equality is due to the fact that 


= A{l_l)(p,g,/Z), 


(E12) 


Thanks to Eqs. (|B3|l and (E17 1 , one gets that the spectral 


decomposition of an Ad^, state with even A can be written 


follows: 


{Cl,C2,C3) (E19) 

= ^[i+c:+(-1)^/2c2 + C3]X|‘1>|X‘I>;I 


2N 

j 

1[1-c:-(-1)^/2^2 + C3]^|<I)-XcI)7I 

j 

^[l+Ci-(-l)^/2^2-C3]Xl'T;)<'I't| 

j 

^[1-c:+(-1)^/2^2-C3]X|'E7X'I'7I. 


Consequently, the spectral decompositions of the Ad^^ states 
represented by the triples ( 5 , 5 X 1 ) and ( 1 , 0 ,/z) are, respec- 
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lively, 

+ jZii+wZiWJi. 

./■ 

and 

^E 21 ) 

j 

+ 5^7(1 


Theorem E.l. For any even N, according to any convex and 
contractive distance, the A1^ state {p, q, 0) is one of the clos¬ 
est {K'^]^^^-separable (resp., M'-separable) states to the 
state (p, q, h). 


Proof. 

Thanks to Lemma E. 1 which holds for any convex and con¬ 
tractive distance and any even N, we just need to prove that for 
any p',q' e [0,1], p' + q' < 1, 


Diip, q, h) , (p, q, 0)) < D((p, q, h) , {p', q',0)). 


In fact 


D((p, q, h) , (p, q, 0)) 


= D 




^ D[Am^(p,q,h),Am^(p',q', 0)) 

< D((p,q,h),{p',q',0)), 


where the first equality is due to Lemma E.2 which holds for 
any contractive distance and any even N, the second equality 
is due to the fact that 


By exploiting the following equalities 

n(|d>})<d>;i) - I't'lX'Pjl, 

Q(|(I)T)<(I)T|) = IT'IX'?}!, 

f^d'i'yx'pyi) = i^yx'pji, 


(E22) 


= (a(E24) 

= (E25) 


with Aji i) representing the LOCC expressed by Eq. (E8i, 
and finally the inequality is due to the contractivity of the dis- 
tance D. 


and the linearity of the channel Q, we immediately get 
Eq. (|E18|l. We then have the inequality 


D({p, q, h) , (p, q, 0)) 




(E23) 


where the final inequality is again due to the contractivity of 
the distance D, whereas the first equality is due to the fact that 




which in turn is due to Eqs. ( |E18| l, ( |E8| l, ( |E5[) and (E4i . By 
putting together the two opposite inequalities ( |Ell| l and ( |E23| l, 
we immediately get the invariance of Eq. ( |E7[ ) for any contrac¬ 
tive distance. 


Now we are ready to find out the analytical expres¬ 
sion of one of the closest {/f^j^j-separable (resp., M'- 
separable) states to any A1^ state vj belonging to the 
{-l,(-l)^''X-l)-corner. 


Now that we know the analytical expression of one of the 
closest and M'-separable states to any AI^ state with 

even N according to any convex and contractive distance, 
we can unveil the general hierarchy of both the and 

M'-inseparable multiparticle entanglement of these AI^ states 
with respect to any geometric entanglement monotone EfX „ 

and respectively. 

Corollary E.l. For every even N and according to any valid 
geometric measure of and M'-inseparable multipar¬ 

ticle entanglement „ and E^,, the following holds: 

^ Up, q, h)) = _ ((p', q, h)) (E26) 

E^fip, q, h)) = EZ,{{p',q, h)) (E27) 

^ Up, q, h)) < _ Up', q, h')), (E28) 

eU(P, q, h)) < EUip',q', h')), (E29) 


for any h < h'. 
Proof. 
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Let us start by proving Eq. ( |E26| l. By using Theorem E.l 
and Lemma [E!21 we obtain 


iip^ q, h)) = D{{p, q, h), (p, q, 0)) 


= D 


1 1 
3’3 


M, 


1 1 
3’3 


,0 


^D((p',q',h),(p\q',0)) 


= E 


((p',q',h)). 


(E30) 


for any p,q,p',q' 6 [0,1 ] p+ q < I, p' +q' < 1 and h 6 [0,1[. 

In order to prove Eq. (E28 1 , let us consider the states 
w = (p,q,h), m' - (p, q, h') such that h < h', and g - - 

g-m' - ip, q, 0), which is one of the closest -s eparable 

states to both m and nr' according to Theorem E. 1 We can 
write m - Am' H- (1 - /l)f, for some d € [0,1]. Now, by using 
the convexity of the distance and Eq. (E26 1 , we get 


E’^^,_^^, iip,q,h)) ^ D{m,g) (E31) 

= DiAm' H- (1 - d)?,?) 

< ADim', ?) H- (1 - A)Dig, g) 

— ADim', g) 

< Dim', g). 


In order to prove Eqs. (E27 1 and ( E29| l, we just remark 
that exactly the same proof holds when substituting 
separability with M'-separability. 


We are now ready to apply the above general results to cal¬ 
culate the geometric multiparticle entanglement E^, ^ (nr) 

and E^Anr) of any Ai^ state nr for particular instances of D. 
As we have just shown, for states in the (-1, (-1)^'^^, -1)- 

corner, £®,|„^(nr) = £",(nr) = D((i, i,/!) ,(i, i,0)) = 

foih), where foW is some monotonically increasing func¬ 
tion of h only, which depends on the chosen distance D. By 
local unitary equivalence, this is true indeed for any Ad^ state 
in any of the four corners if we define the generalised h to be 
hju - \iYA’j=i k;l - 1) ■ Table I in the main text shows foihm) 
for the relative entropy, trace, infidelity, squared Bures, and 
squared Hellinger distance. 

Here we show the derivation of the given expressions for 
foih). Since the two Ad^ states ( 5 , 5 , /r) and Q, j, o) are diag¬ 
onal in the same basis, we have that their distance reduces to 
the corresponding classical distance between the probability 
distributions formed by their eigenvalues, denoted by f*/, and 
Pq respectively. We recall that the classical relative entropy, 
trace, infidelity, squared Bures, and squared Hellinger dis¬ 
tance between two probability distributions P = jp, ) and Q — 
[qi) are given by, respectively Dre(P, Q) = TjiPi\og 2 iPilq,), 
DtAP, Q) = i:,- \Pi - qi\ /2, D^iP, 2) = 1 - FiP, Q), and 
Db(P,0 = DniP,Q) = 2(1 - ^mp;W), where FiP,Q) = 


(Z/ -\/Piqi) is the classical fidelity. Consequently, by using 
Eq. ( |E20| i to get both P/, and Pq, we obtain the desired expres¬ 
sions for faih). 


2. Odd N case 


Let us now turn our attention to the evaluation of the geo¬ 
metric multiparticle entanglement Pf* - im) and E^Am) of 

an Al^ state m in the case of odd N. Unlike the even N 
case, where all the convex and contractive distances D con¬ 
cur on what is one of the closest separable states to an Al^ 
state, there is unfortunately no such agreement in the odd N 
case. In the following we will focus in particular on the trace 
distance-based geometric measures of multiparticle entangle¬ 
ment i'^) E^; im) of any AI^ state with odd N, 

when considering as usual a partition such that K'^ 

is odd for more than one value of a and a number of parties 
M' > \NI2\ respectively. 

We know that the trace distance-based multiparticle entan¬ 
glement of an Ad^ state m with odd N is the minimal dis¬ 
tance from m to the unit octahedron 0\. Due to convexity 
of the trace distance and of the unit octahedron 0\, we get 
that one of the closest -separable (resp., M'-separable) 

Ad^ states to an entangled Ad^ state m must belong necessar¬ 
ily to the boundary of the octahedron, i.e. either to one of its 
faces or edges. We can easily see that the trace distance be¬ 
tween two arbitrary odd N states nri = |cj^\ \ j and 

m 2 - is nothing but one half of the Euclidean 

distance between their representing triples, i.e. 


Djrimi,m2) “ 2 ^ 




(E32) 


This is proven as follows. In order to evaluate the trace dis¬ 
tance between any two Ad^ states with odd N, we just need to 
calculate the eigenvalues of their difference, because 

DTi-(riTi,'nT2) = ^Tr(|uTi - m2\) ^ 


where {T,) are just the eigenvalues of mi - m 2 - Since 
'^i - '^2 - 2 V Tii d-icrf^, with dj - - cf \ one can eas¬ 

ily see that its eigenvectors are exactly the ones expressed 


in Eq. (B5i, with di - rsinflcos^, d 2 - rsin0sin0 and 


do, - rcosO while its eigenvalues are given by either or 

-^r, with r = .^d^ ‘^2 putting these eigenvalues 

into Eq. (E33 1 one immediately gets Eq. (E32 1 . 

An immediate consequence of Eq. (E32 1 is that the clos¬ 
est -separable (resp., M'-separable) Ad^ state = 

[ii, S 2 , ^ 3 ) to an entangled Ad^ state m - {ci,C 2 ,C 3 ) is just its 
Euclidean orthogonal projection onto the boundary of the unit 
octahedron Oi- We can thus distinguish between the follow¬ 
ing two cases; 
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1. the Euclidean projection of {ci,C 2 ,C 3 ) onto the bound¬ 
ary of the unit octahedron 0\ falls onto one of its faces. 
This case happens if, and only if, 0 < sign(c,)i,' < 1 for 
any i, where {i/ = sign(c,)(l - |ci| - |c 2 | - ksl -H 3|c/|)/3) 
is exactly the triple representing such Euclidean projec¬ 
tion; 

2. the Euclidean projection of {ci,C 2 ,C 3 ) onto the bound¬ 
ary of the unit octahedron Oi falls onto one of its edges. 
This case happens when any of the conditions listed in 
case Edo not hold. Moreover, the triple {ii, S 2 , 53 ) rep¬ 
resenting such Euclidean projection is given by = 0 
and Si = sign(c,)(l k^j -i- 2 |c,j)/ 2 , where k is set by 

fk = min{/i,/ 2 ,/ 3 ) with/- = + (l - 

This provides the explicit expression reported in Eq. (3) in the 
main text. 


Appendix F: Genuine geometric muitiparticie entangiement 
of the GHZ-diagonai states 

In this appendix we will adopt our approach to evaluate ex¬ 
actly the geometric genuine multiparticle entanglement of 
any A^-qubit GHZ-diagonal state ^ with respect to any con¬ 
tractive and jointly convex distance D. 

Recall that any A^-qubit state g can be transformed via a 
single-qubit LOCC E into a GHZ-diagonal state gcnz = r(p) 
with eigenvalues given by pf = (J3f\g\/3f} IfTTI . a procedure re¬ 
ferred to as GHZ-diagonalisation of g in the main text. The en¬ 
tanglement quantihcation is then based on the following two 
arguments. 

Eirst, we have that one of the closest 2-separable states to 
a GHZ-diagonal state is itself GHZ-diagonal. Indeed, for any 
GHZ-diagonal state ^ and any 2-separable state g, we have 
that 


D(^, gcHz) = Dim, r(?)) < Di^, g), (El) 


Now, let us consider an arbitrary GHZ-diagonal state ^ and 
rearrange its GHZ eigenstates in such a way that the 

corresponding eigenvalues are in non-increasing order. 

It is well known that ^ is 2-separable if, and only if, <1/2 
0. Eor pj > 1/2, we will show that one of the closest 

2-separable GHZ-diagonal states g^ has eigenvalues {pf}^^^ 
such that pf = 112, with {pf}^!^ corresponding again to the 
ordering of GHZ eigenstates set by This result fur¬ 


ther simplihes the optimisation in Eq. E2 


Consider any 2-separable GHZ-diagonal state g, it holds 
that there will always be a 2-separable GHZ-diagonal state g' 
with eigenvalues {p[ and pj = 1/2 such that g' - + 

(1 - /!)§• for some A e [0,1]. Now, for any convex distance, 
the following holds 


D(^,g') (E3) 

^D(^,A^ + il-A)g) 

<ADi^,^) + (l-A)D(^,g) 

^(l-A)Di^,g) 

<Di^,g). 

This inequality immediately implies that one of the closest 
2-separable GHZ-diagonal states g^ to a 2-inseparable GHZ- 
diagonal state ^ is of the form g', which is formalised as a 
corollary below. 

Corollary F.l. For any convex and contractive distance D 
and any fully inseparable GHZ-diagonal state f, whose GHZ 
eigenstates are arranged such that the correspond¬ 

ing eigenvalues {pf)?”[ are in non-increasing order, one of the 
closest 2-separable states g^ to ^ is itself a GHZ-diagonal 
state with eigenvalues {pf}f^ such that pf — 1 / 2 , with 
{pf}fli corresponding again to the ordering of GHZ eigen- 
states set by f. 


where in the first equality we use the invariance of any GHZ- 
diagonal state through E and that r(^) = gcnz is the GHZ- 
diagonalisation of g, and in the inequality we use the contrac- 
tivity of the distance through any completely positive trace¬ 
preserving channel. Moreover, the GHZ-diagonalisation ^ghz 
of any 2-separable state ^ is a 2-separable GHZ-diagonal state 
since E is a single-qubit LOCC. 

Therefore, the set Sf of 2-separable GHZ-diagonal states 
turns out to be the relevant one in order to compute exactly 
any distance-based measure of genuine multiparticle entan¬ 
glement of a GHZ-diagonal state f, thus dramatically simpli¬ 
fying the ensuing optimisation as follows: 

Ef(f)= inf D(f,g)^ inf D{f,gcnz). (F2) 

fGHze^ 


We can now apply this Corollary to calculate the geomet¬ 
ric genuine multiparticle entanglement Ef(f) of any GHZ- 
diagonal state f for particular instances of D. Since the closest 
2-separable state g^ to a GHZ-diagonal state f is also a GHZ- 
diagonal state, they are diagonal in the same basis and their 
distance reduces to the corresponding classical distance be¬ 
tween the probability distributions formed by their eigenval¬ 
ues, denoted by and respectively. By using the expres¬ 
sions given earlier in Appendix]^ of the classical relative en¬ 
tropy, trace, inhdelity, squared Bures, and squared Hellinger 
distance between two probability distributions P^ and P^^, 
and minimising it over all probability distributions P^^ such 
that pf - 1 / 2 , one easily obtains the desired expressions for 
Efif) expressed in the main text. 




